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Abstract 

The effect of the Coulomb-interaction on persistent currents in disordered 
mesoscopic metal rings threaded by a magnetic flux cj) is studied numerically. 
We use the simplest form of self- consistent Hartree theory, where the spatial 
variations of the self-consistent Hartree potential are ignored. In this approxi- 
mation the self-consistent Hartree energies are simply obtained by diagonaliz- 
ing the non-interacting system via the Lanczos method and then calculating 
the (disorder-dependent) particle number on the ring self-consistently. In 
the diffusive regime we find that the variance of the total particle number is 
strongly reduced, in agreement with the prediction of the random-phase ap- 
proximation. On the other hand, the variance of the number of energy levels 
in a small interval below the Fermi energy is not affected by the Coulomb 
interaction. We argue that this implies that the experimentally observed en- 
hancement of the persistent current is due to long-range Coulomb interactions. 

PACS numbers: 73.50.Bk, 72.10.Bg, 72.15.Rn 
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I. INTRODUCTION 



The usual diagrammatic approach [|IJ] to disordered metals is based on the existence of 
two small parameters: the smallness of (kpi)" 1 (where kp is the Fermi wave-vector and I 
is the elastic mean free path) justifies the perturbative treatment of the disorder potential, 
and the smallness of the parameter njkp at high densities (where k is the Thomas- Fermi 
screening wave- vector) justifies the so-called random-phase approximation (RPA), which re- 
sums the most divergent diagrams in the expansion of the Coulomb potential. This leads to 
the screening of the long-range tail of the Coulomb potential, so that the electron-electron 
interaction can be taken into account via an effective short-range interaction of the Hubbard- 
type. 

In Ref. [0 this standard diagrammatic approach was used to calculate the effect of 
electron-electron interactions on the average persistent current in mesoscopic metal rings 
threaded by a magnetic flux cf) Surprisingly, the average current was found to be 

almost two orders of magnitude smaller than the experimental result by Levy et al. ||. This 
experiment and subsequent experiments by other groups ||[7j have stimulated many recent 
theoretical works. Nevertheless, till now a generally accepted theoretical explanation of 
the experimental data |||| has not been found, perhaps because the standard perturbative 
approach is not applicable in the persistent current problem. This has several reasons: 
First of all, the experiments [|5]|| are performed at constant particle number N, so that the 
persistent current should be calculated from the flux-derivative of the canonical free energy 
F(N,<p), 

where ip = 4>/4>q is the flux measured in units of the flux quantum 0o — hc/e. On the other 
hand, with the standard machinery of many-body theory one calculates the current /(//, ip) 
at constant chemical potential /i, which is related to the grand canonical potential fi(/z, p) 
via 

edQ(n,(p) 

1{ ^ ) = -h^p— ■ (2) 
Several authors pHTlH have noticed that for non-interacting electrons there exist striking 



differences between the average currents I(N, p) and /(/i, <p) (the overline denotes averaging 
over the disorder). In particular, for non-interacting electrons in the diffusive regime I(fi, p) 
is exponentially small, whereas I(N, p) is not. The question whether these qualitative 
differences persist even in the presence of electron-electron interactions has not yet been 
investigated. 

We would like to point out two more possible reasons for the failure of conventional per- 
turbative many-body techniques in the persistent current problem. The first is the existence 
of two different purely geometric length scales in mesoscopic rings, namely the circumference 
L, and the characteristic thickness L±_ of the rings. Thus, diagrams which are usually ne- 
glected because they involve higher orders in {kpt)~ l may become important because they 
are geometrically enhanced by additional powers of Lj 'L± ^> 1. This has first been pointed 
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out by Beal-Monod and Montambaux Finally, it is questionable whether the conven- 
tional RPA is valid for persistent currents. Note that the flux- dependent parts of F(N, ip) 
and ip) are much smaller than the flux- independent parts. The RPA has been developed 
to re-sum the leading (flux-independent) terms in the expansion of f2(/z, <p) in powers of the 
Coulomb interaction However, it is not at all clear whether such a resummation of a 
formally divergent series remains valid for the sub-leading flux- dependent part of Q([i,<p). 
We shall come back to this point in Sec.|TTJ. 



Recently one of us fb4]J15| has proposed that the flux-dependent part of the average 
Hartree energy in a mesoscopic ring is only weakly screened, and showed that in this way 
the large magnitude of the experimentally observed currents |j| can be easily explained. 
The arguments given in Refs. jl4|Jl5| were simple but powerful, and made use of the exact 
Hohenberg-Kohn theorem of density-functional theory ||16|| . Although Ref. |TJ]] has been 



criticized [|17yi8| , all arguments put forward against Ref. |L4j where essentially based on the 
RPA. The possibility that the RPA might not be a good approximation for the calculation 
of the persistent current [19|]2(| was ignored. Because even in the absence of disorder it is 



very difficult to calculate the dielectric function beyond the RPA, numerical methods seem 
to be the only controlled way to settle the issue. So far all numerical studies of two- and 
three-dimensional systems in the diffusive regime have found that the long-range part of 
the Coulomb interaction indeed strongly enhances the average persistent current |21|-p4 



Among the numerical methods, the exact diagonalization of the many-body Hamiltonian 



21| , |22[| is free of any approximations. Unfortunately, it is very difficult to study systems 
with more than approximately 10 electrons with this method. This is certainly too small 
to address the screening problem. Calculations based on the self-consistent Hartree-Fock 
approximation have been pushed to electron numbers of the order of 10 2 . Although these 
calculations do not produce the exact solution of the many-body system, they are non- 
perturbative by virtue of the self-consistency condition, and are therefore more reliable than 
any low-order diagrammatic calculation. 

Guided by the evidence that the enhancement of the persistent current is closely related 
to the physics of screening, we shall in this work examine the screening problem numerically. 



As a starting point, we follow Ref. Jy] and use the approximate mapping from the average 
canonical current onto an effective grand canonical average, 

WM-TWM^^^^M , (3) 



where I(fi*,ip) is the average grand canonical current at chemical potential fi*, and 



kM J^eS) (4) 



is the average compressibility. £ 2 (/i*,y?) is the variance of the particle number N([i*,(p) in 
a grand canonical ensemble 
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Y, 2 (jS,<p) = N*(ji;<p)- [N(fi*,<p)\~ . (5) 

The value of fi* should be chosen such that the disorder- and flux-averaged particle number 
in the corresponding grand canonical ensemble agrees with the given particle number N in 
the original canonical ensemble, i.e. 
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/ dipN(u,*, Lp) = N 
Jo 



(6) 



Note that for non- interacting electrons in the diffusive regime <p) is exponentially small 



25 1 . The leading interaction contribution to I(fM*,(p) has been calculated by Ambegaokar 



and Eckern ||. In this work we would like to focus on the effect of electron-electron interac- 
tions on the difference between the canonical and the grand canonical current, i.e. the term 
on the right-hand side of Eq.|3]. Recently Berkovits and Avishai |22| studied this term via 
exact diagonalizations of small systems of electrons interacting with Coulomb forces. How- 
ever, with exact diagonalizations it is impossible to reach system sizes where the condition 
£ <C L <C £ is realized. Here t is the elastic mean free path, L is the circumference of the ring, 
and £ is the localization length. This condition defines the diffusive regime, which is relevant 
for the experiments of Refs. In this work we shall therefore take a different approach, 

and evaluate the right-hand side of Eq|| within the simplest possible approximation which 
still contains the physics of screening, namely the simplified capacitance model |15|J20|J26|1 . In 
this model the spatial dependence of the self-consistent Hartree potential is ignored, so that 
it simply renormalizes the chemical potential. Nevertheless, diagrammatically this model 
still contains the usual infinite series of bubble diagrams which lead to the RPA-picture of 
screening. But the model contains infinitely many other non-RPA diagrams, so that the 
numerical solution of the self-consistent Hartree theory for this model allows us to check the 
accuracy of the RPA. Note that the simplified capacitance model is easily solved numerically 
by diagonalizing the non-interacting system with the help of the Lanczos method, and then 
imposing a simple self-consistency loop for the particle number. The most difficult step in 
this procedure is the exact diagonalization, so that in this way we can reach the same system 
sizes as without interactions. 



II. THE SIMPLEST FORM OF SELF-CONSISTENT HARTREE THEORY 

To calculate the effect of electron-electron interactions on the persistent current, let us 
use the self-consistent Hartree approximation to take the Coulomb interactions between all 
charges (i.e. electrons and positively charged ions) into account. Numerically it is convenient 
to use a lattice model. The Hartree energies E a ((p) and wave-functions ^ Q (r) are obtained 
by solving the Schrodinger equation 



E 

fi=x,y,z 



[t^ a (r + a M ) + t*^ a (r - a„)] + [U(r) + U H (r)} ^ Q (r) = E a { V )tfj a {r) , (7) 



where r labels the sites of a three-dimensional cubic lattice with primitive vectors & x , & y , 
a 2 . To model the Aharonov-Bohm flux through the ring, we call the azimuthal direction 

the x-direction and choose t x = te N * , t y = t z = t. Here N x is the number of lattice sites in 
the x-direction. For convenience, from now on all energies will be measured in units of the 
hopping energy t. As usual, disorder is introduced via random potentials U(r), which are 
assumed to be independent random variables with zero average and uniform distribution 
in the interval [— w/2,w/2]. The self-consistency is imposed by requiring that the Hartree 
potential Uh(t) satisfies 
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e 2 



M^E^KO-n+M] , (8) 

r' 



r 



with n(r) = X^'IV^ 1 ")! 2 - The prime indicates that the a-sum is in a canonical ensemble 
over the lowest N solutions of Eq.^, and in a grand canonical ensemble over all a with 
E a < //, where \x is the chemical potential. The number density n + (r) of the positive 
background charge is assumed to be fixed. Given the self-consistent solutions of Eqs.[7| and 
H the canonical persistent current can be calculated from Eq.|l] with [16 



N x 

F(N,<p) = ^E a (<p)--^n(T)U s (T) , (9) 

a=l Z r 

and the grand canonical current is obtained from Eq.[| with Q(/i,ip) = F(N(n,ip),ip) — 
fj>N(n, <f), where the function N(p, ip) is now explicitly given by 

N0i,ip) = -E a (ip)) . (10) 

a 

In order to obtain the average current, we should solve Eq.f7| self-consistently for many 
realizations of disorder and then average the result. Of course, this is a very difficult problem, 
which is beyond the scope of this work. Instead, we shall use Eq.||| to calculate the difference 
between the canonical- and the grand canonical current approximately. 

To reduce the numerical work even further without loosing the physics of screening, let 
us replace the self-consistent Hartree potential by its spatial average, Uh(t) — > ^ J2r Un{ r ), 
where N L = Y, r * s the number of lattice sites fl5p0p6H . In this approximation Uh is simply 
given by 



1 H 



<iN - N+) 1 1 Ei-^-77 • (11) 



Co Co N L (N L - 1) ^ |r 

r^r' 



where iV = J2 r n ( r ) is the total number of electrons on the ring, and iV + = J2 r n +( r ) 
is the total number of positively charged ions. Note that Co is nothing but the classical 



capacitance of the ring. In Fourier space the approximation |TT| means that we neglect the 
momentum transfer through all Coulomb lines. Because for a thin ring (with Lj_ L) the 
long-wavelength Fourier components of the Coulomb potential depend only logarithmically 
on the wave-vector, we expect that in this case the approximation |TT] is sufficient for a 
qualitative estimate of the importance of Coulomb interactions. Of course, in a canonical 



ensemble Eq.[TT]is a trivial flux-independent constant, so that the interaction-contribution to 



the persistent current vanishes. For this reason only the non-zero Fourier components of the 



Coulomb-potential were retained in Ref. At long wavelengths these can be expressed in 
terms of generalized capacitances C/., describing long- wavelength charge fluctuations which 
do not change the total electron number. On the other hand, in a grand canonical ensemble 
the particle number iV(/i, ip) is a very complicated flux- dependent random variable, which 



has to be determined self-consistently by solving EqJlOL By making the approximation 11 



on the right-hand side of Eq.|3], we can estimate the interaction contribution to the canonical 
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persistent current in a way that is very well suited for numerical calculations. For the 
simplified capacitance model Eq.10 reads 



Co ' 



e 2 iV 4 



(12) 



where the energies e a ((p) are the solutions of Eq.[7| for Uh = 0. Eq.[L2] can be viewed as the 
simplest form of self-consistent Hartree theory. 

To appreciate the non-perturbative nature of Eq.|12|, let us calculate the average com- 
pressibility k(/i, ip) defined in Eq||. Therefore we simply differentiate both sides of Eq.|l2| 
with respect to // (taking into account that dN + /dfi = 0, because the positive background 
charge is fixed), solve for dN/dfi, and finally average. This yields 



1 + TTp(P^) 
^0 



(13) 



where the self-consistent density of states is given by 



P0,^>) = ^6{e-e a {<p) - 



e 2 (jV(/z^)-iV 4 
Co 



(14) 



Ignoring random fluctuations of the density of states, we may factorize the average, 



p(p,<f) 



1 + 7^P(/W) 



(15) 



In the diffusive regime it is also reasonable to neglect interaction contributions to the average 



density of states | 26fl , so that 



P(P, <p) Po(p, <p) = ^2S{p- e a (ip)) 



A 



(16) 



where A is the average level spacing at the Fermi energy in the absence of interactions. 
Hence we obtain 



>RPA 



A" 



1 



'RPA 



1 + ^ 



(17) 



Note that A -1 is the average compressibility without interactions. Because e 2 / (CoA) ^> 1 
in the experimentally relevant parameter regime, the RPA renormalization factor -Zrpa 
is small compared with unity. Hence, the interactions lead to a drastic reduction of the 
average compressibility. A similar result has also been obtained by Vignale within a density 
functional approach |S7| . A strong reduction of the compressibility has also been observed 
in the exact diagonalization study by Berkovits and Avishai [p2]| . The origin for this effect 
is intuitively obvious, because physically k _1 is the energy associated with the addition of 
an electron to the system. In the absence of interactions this is simply the average level 
spacing, but with Coulomb interactions one has to pay the charging energy e 2 /C*o for adding 
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an electron. Hence, the inverse compressibility is pushed from A to the much larger energy 
e 2 /Cq. We would like to stress that in the conventional diagrammatic approach to the many- 
body problem for infinite translationally invariant systems the Hartree diagram responsible 
for this renormalization is usually ignored, because it is assumed to be canceled exactly by 
the positive background charge p5[| . Substituting Eq.[I7| into EqfJ, we obtain 

~t4| e V,,) , (18) 

where the second line is valid for e 2 /Co ^> A. The crucial question is now whether the 
flux-dependent part of the variance £ 2 (/i*,<£>) is reduced by the interactions or not. Note 
that naive application of the standard RPA |[7|,iZl leads to the replacement £ 2 (/i*,y?) — > 



Z^ PA EQ(/i*, if), where T,l(fi*,ip) is the variance in the absence of interactions (see Eq.|20| 
below). However, as already mentioned, the flux- dependent pa.it of S 2 (/i*, ip) is much smaller 
than the flux-independent part, and it is not at all clear whether it is screened in exactly the 
same way as the flux- independent part. Because even without disorder it is very difficult to 



calculate corrections to the RPA in a systematic way ||13|| , we shall in this work study this 
problem numerically. 



III. NUMERICAL RESULTS 



We have calculated the eigenvalues E a {(p) = e a (ip)- s r^(N(fj 1 ,ip)—N + ) of Eq.[7] numerically 
for fixed realizations of the random potential U(r) on finite lattices. The eigenvalues e a (ip) 



without the Hartree term were calculated with the Lanczos method To obtain the 

energy shift due to the Hartree potential, the self-consistency equation [T2| for the particle 
number for a given realization of the random potential was then solved via a simple Newton 
procedure. It turns out that at zero temperature this equation does not always have a 
unique solution, because by definition N must be an integer. This problem is easily avoided 
by working at finite temperature T, and then extrapolating for T — > 0. Because for T > 
the step function in Eq.|l^ is smoothed out into the Fermi function, there exists always a 
unique solution of Eq.[L2|. This is easy to see from the fact that the left-hand side of Eq.|l2| 
is a monotonically increasing function of the particle number, while the right-hand side is 
monotonically decreasing. 

Most theoretical calculations in disordered systems are based on the assumption that in 
the energy window of interest the energy-dependence of the average density of states can be 
ignored. On the other hand, in our finite system the average density of states is non-zero only 
in a finite interval, and exhibits a broad maximum at zero energy, as shown in Fig.|I| As we 
have recently pointed out for non-interacting electrons the universal weak localization 



effects related to the Cooperon pole are completely washed out if the energy- dependence of 
the average density of states becomes significant in the energy interval under consideration. 
Hence, in order to compare our numerical calculations with theoretical predictions for models 
with constant density of states, we shall retain only the central part of the spectrum, and 
discard all energies outside an interval [/i — E , (/]. From Fig.[I| it is obvious that for our set 
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of parameters Eq = 2 is a good choice. A similar procedure has also been adopted in Ref. 



32fl . Note that in this way we formally replace the density of states of our original tight 
binding model by the density of states with sharp edges shown in Fig.|I|. The total particle 
number is still defined as in Eq.^, except that now the a-sum is over energies satisfying 
e a G — E ,fj]. The great advantage of this construction is that the average density of 
states is practically constant in the entire regime where it is non-zero, so that non-universal 
effects related to the energy-dependence of the density of states should be clearly visible at 
the band edges. 

The drastic effect of the Coulomb interaction on the statistics of the number of electrons 
in the system is clearly seen in Fig.|2], were we plot the probability distribution P{N) of the 
particle number for different values of e 2 /Cq. The narrowing of the distribution for increasing 
interaction parameter can be viewed as a screening effect. The usual RPA prediction for the 
second moment £ 2 (/x, </?) = (iV(/i, <£>) — N + ) 2 of the distribution P(N) around N + is easily 
obtained from Eq.|l2| by expanding the right-hand side to first order in N — N + , 

N-N + = N -N + -Z r p { l i,<p)(N-N + )+O{(N-N + ) 2 ) , (19) 

where N = J2 a ©(A* ~ e a) is the particle number in the absence of interactions, and po(p, <p) 
is the non- interacting density of states, see Eq.[T^. Solving Eq.|TU| for iV — N+, squaring, 
averaging, and ignoring random fluctuations of the density of states (see Eq.[L5|), we obtain 
for the second moment 

£W) = Z 2 PA £ 2 (/i,^) , (20) 



where Zrpa is given in Eq.[L7|, and £q(/z,<^) = (N (p,ip) —N + ) 2 . The obvious question is 
now whether Eq.BOl is accurate or not. Note that this expression is based on the expansion 



19| , which is a priori uncontrolled, because it is not clear whether the effect of the higher 
order terms that have been ignored is important or not. Furthermore, it is by no means 
clear whether the flux-dependent part of £ 2 (/U, <p) is modified in exactly the same way by 
the interaction as the dominant flux-independent part. To investigate this point, let us 
expand £ 2 (/i, <p) in a Fourier series. Because N(p, <p) is an even periodic function of if with 
a fundamental period of unity, we may expand 

p(0) oo 

EW) = — + ]TpWcos(27my,) . (21) 

A n=l 

Denoting by P the corresponding Fourier components in the absence of interactions, the 
RPA predicts that pM/P (ri) = Z£ PA for all n. As shown in Fig.|3], for the zeroth and the first 
two even Fourier components (n = 2m = 0,2,4) this prediction is in excellent quantitative 
agreement with our numerical self-consistent solution of Eq.|12l The odd components (n = 



2m + 1) are numerically found to be much smaller than the even ones |pl |, but are screened 
in a similar fashion. We conclude that for our simplified capacitance model the RPA is an 
excellent approximation, so that the interaction actually reduces the current in Eq.|l8| by a 
factor of Zrpa <C 1. In this case the difference between the canonical and the grand canonical 
current is completely negligible compared with the leading contribution to the average grand 
canonical current I(p*,(p) that has been considered by Ambegaokar and Eckern [H. 
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However, the screening behavior is completely changed if we exclude the energies in the 
vicinity of the band edges. Such a procedure can be motivated physically as follows: It is 
well-known that in an interacting many-body system propagating quasi-particles exist only 
for energies sufficiently close to the Fermi energy. Moreover, the flux-dependence of physical 
observables should be due to quasi-particles that can coherently propagate around the ring, 
and in this way probe the sensitivity to twists in the boundary conditions. Hence, we expect 
that in a realistic interacting many-body system <9£ 2 (/i, if)/dif is essentially determined 
by weakly damped states with energies in a small interval of width E in below the Fermi 
energy, Here E- in is some inelastic cutoff energy that is small compared with the Fermi 
energy. In other words, inelastic processes (which are completely ignored in the Hartree 
and the Hartree-Fock approximation) will naturally restrict the states contributing to the 
flux- dependent part of S 2 (-E, f) to those with energies in the interval \p — E- m , /i]. The theory 

is based on the existence of such a cutoff energy. 



developed in Refs. p^jll 

To study the contributions from the part of the spectrum close to the Fermi energy to 
T, 2 (p,ip), let us therefore introduce the auxiliary quantity 



a 2 (E,p,<f) 



/ del de'K 2 (e,e',cp) 

Ja-E Ju—E 



(22) 



K 2 (e, e', <f) = p(e, <f)p(e', <f) - p(e, if) p(e', if) , (23) 

where the self-consistent density of states p(e,<p) is given in Eq.[l4]. Note that a 2 (E,p,<f) 
is the variance of the number of energy levels in an interval of width E below the Fermi 
energy, and that by construction £ 2 (/i, <f) = lim^^oo <? 2 (E, p, if). In fact, given our truncated 
density of states shown in Fig.|I], we have £ 2 (/i, <p) = cr 2 (E, p, if) for all E > E . It is again 
convenient to expand cr 2 (E, p,if) in a Fourier series, 



a (E,n,tp) 



P 



+ Y J P {n \ E ) cos(27mv?) 



(24) 



n=l 



In the absence of interactions and for energies E small compared with H/r (where r is the 
elastic lifetime) an approximate expression for the Fourier coefficients Pq\e) (where the 
subscript indicates that interactions are neglected) can be obtained from the Feynman 
diagrams with two Cooperons and two Diffusons given by Altshuler and Shklovskii ||33|| . In 
this approximation one obtains p l Q Tn+1 \E) = 0, and the even Fourier components are given 



by 3| 



P ( o m \E) 



7r 2 m 



exp 



-mt 1 ' 2 



cxp 



X COS 



m 
~V2 
m 



E 2 



„\ 1/2' 

r 



_\ 1/2 

r 



(25) 



The flux average is 



P$\E) 



7T- 

oo 



Erf m) (£) 



e 2 + r + r 



fl/2 



(26) 



m=l 
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Here E = E/E c , and T = T/E c , where E c is the Thouless energy and T is some cutoff 
energy that has been introduced into the non- interacting theory by hand ||. For non- 
interacting electrons on a quasi-one dimensional ring (where diffusion is only possible along 
the circumference) we have recently confirmed Eqs.^5| and |26] numerically Here we 

would like to investigate whether in the presence of long-range Coulomb interactions these 
expressions are strongly reduced by screening effects. Note that according to Altland and 
Gefen (18j there should be no qualitative differences between the screening of fluctuations 
of the spectral density and the total density; if this would be correct, then the variance 
a 2 (E, fi,ip) of the number of energy levels in an interval E below the Fermi energy should 
be screened in precisely the same way as the variance £ 2 (/i, <fi) of the total particle number 
(which is strongly affected by screening, see Fig.|||). We are now in the position to settle 
this controversy numerically within the simplified capacitance model. We would like to 
emphasize that, in spite of its simplicity, this model contains the physics of screening. 

In Fig.[| we show our numerical results for the zeroth Fourier component p^ (E) as func- 
tion of E for different values of e 2 /(C*oA). Corresponding results for the first even Fourier 
component p^ 2 \E) are shown in Fig.|5]. Evidently, in the entire energy interval where the 
average density of states is approximately constant, the flux average of <r 2 (E, \x, <p) and the 
dominant flux-dependent contribution are only weakly affected by the Coulomb interac- 
tion. Moreover, in Fig.|| we also show the Altshuler-Shklovskii prediction [25| for p^ 2 \E). 
The agreement with the numerical results for the simplified capacitance model clearly 
demonstrates that for small energies it is indeed allowed to ignore screening corrections 
to cr 2 (E, /i, if). Thus, our numerical results support the arguments put forward by one of us 
in Refs. |ri^Jl5|JT9|j20[| , and strongly disagree with the statement of Altland and Gefen [18| 
that the fluctuations of the spectral density is screened in the same way as the total density. 
From Figs.f| and [5] it is also obvious that the RPA screening of the fluctuation of the total 
particle number is recovered as soon as the interval E includes the band edge Eq. The sharp 
drop for p(°\E) and p^ 2 \E) for E m Eq is in agreement with the RPA result for the Fourier 
components of the variance of the total particle number shown in Fig.^. Of course, to 
obtain the physical persistent current, fluctuations on all energy scales should be included, 
so that the simplified capacitance model does not show any enhancement of the average 
persistent current due to long-range Coulomb interactions. However, as explained in Refs. 



|T5| , |3T|| , we expect that in a realistic many body system the damping of quasi-particles which 
are not close to the Fermi surface will eliminate the contribution from high-energy states to 
the flux- dependent part of £ 2 (/i,y?). In other words, for a realistic interacting many-body 
system we expect 

— £ 2 (/i,<^) ^ — a 2 (£,/i,<^) , for all E > E- m , (27) 



where E in is the unknown inelastic cutoff energy mentioned above. Note that in Ref. fTT 



Altshuler, Gefen, and Imry implicitly seem to make a similar assumption 34 



IV. CONCLUSIONS 

In this work we have presented numerical results for the simplified capacitance model, 
which describes the effect of long-range Coulomb interaction on the average canonical per- 
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sistent current. The model can be considered as the simplest form of self-consistent Hartree 
theory, which still contains the non-trivial physics of screening. We have confirmed one of 
the fundamental assumptions of Ref. [ Hf , namely that the fluctuations of the spectral den- 



sity for sufficiently small energies are not screened. Of course, within the approximations 
inherent in our model all eigen-energies (including those far away from the Fermi energy) 
correspond to quasi-particles with infinite lifetime, so that the physical persistent current is 
eventually screened once the energies at the band edges are taken into account. However, 
in physically more realistic models we expect that the damping of the quasi-particles will 
automatically eliminate non-universal effects related to the band edges. Given the fact that 



the low-energy fluctuations are not screened, the universal current proposed in Ref. 
directly follows. 

The problem of verifying by explicit calculation that the damping of quasi-particles 
restricts the flux-dependent part of the spectrum to a small interval of energies close to the 
Fermi energy remains open. We believe that this problem is closely related to the rather 
obscure inelastic cutoff T, which is usually introduced phenomenologically into the Cooperon 
and Diffuson propagators of the non-interacting theory. For a recent microscopic calculation 
of r see Ref. |35| . More generally, the inelastic cutoff should be momentum- and frequency- 
dependent, r(q, uS). It is tempting to speculate that the energy scale where r(0, u) begins 
to deviate significantly from its zero-energy limit T(0, 0) is proportional to the cutoff energy 
E- m mentioned above (see Eq 



ACKNOWLEDGEMENTS 



We are grateful to J. A. Holm for providing us with his Master-Slave Implementation for 
Parallel Virtual Machines [30|, which opened the way for an efficient numerical simulation 
on a work-station cluster. We would also like to thank K. Schonhammer for discussions. 
This work was financially supported by the Deutsche Forschungsgemeinschaft (SFB 345). 



11 



REFERENCES 



[1] See the articles by B. L. Altshuler and A. G. Aronov, or by H. Fukuyama, in Electron- 
Electron Interactions in Disordered Systems, eds. A. L. Efros and M. Pollak, (North 
Holland, Amsterdam, 1985), and the review by P. A. Lee and T. V. Ramakrishnan, 
Rev. Mod. Phys. 57, 287 (1985); 
[2] V. Ambegaokar and U. Eckern, Phys. Rev. Lett. 65, 381 (1990); ibid. 67, 3192 (1991); 
U. Eckern, Z. Phys. B 82, 393 (1991). 
F. Hund, Ann. Phys. (Leipzig) 32, 102 (1938). 

M. Biittiker, Y. Imry, and R. Laudauer, Phys. Lett. 96 A, 365 (1983). 
L. P. Levy et al, Phys. Rev. Lett. 64, 2074 (1990). 
V. Chandrasekhar et al, Phys. Rev. Lett. 67, 3578 (1991). 
D. Mailly, C. Chapelier and A. Benoit, Phys. Rev. Lett. 70, 2020 (1993). 
H. Bouchiat and G. Montambaux, J. Phys. (Paris) 50, 2695 (1989). 

A. Schmid, Phys. Rev. Lett. 66, 80 (1991). 

F. von Oppen and E. K. Riedel, Phys. Rev. Lett. 66, 84 (1991). 

B. L. Altshuler, Y. Gefen, and Y. Imry, Phys. Rev. Lett. 66, 88 (1991). 
M. T. Beal-Monod and G. Montambaux, Phys. Rev. B 46, 7182 (1992). 
See, for example, R. D. Mattuck, A Guide to Feynman Diagrams in the Many-Body 
Problem, 2nd ed. Chap. 10 (Dover, New York, 1992; 1st ed. by McGraw-Hill, New York, 
1967). 

P. Kopietz, Phys. Rev. Lett. 70, 3123 (1993); ibid. 71, 306 (1993) (E). 
P. Kopietz, Int. J. Mod. Phys. B 8, 2593 (1994). 
P. Hohenberg and W. Kohn, Phys. Rev. 136, B864, (1964). 

G. Vignale, Phys. Rev. Lett. 72, 433 (1994). 
A. Altland and Y. Gefen, Phys. Rev. Lett. 72, 2973 (1994). 
P. Kopietz, Phys. Rev. Lett. 72, 434 (1994). 
P. Kopietz, Phys. Rev. Lett. 72, 2974 (1994). 
R. Berkovits and Y. Avishai, Europhys. Lett. 29, 475 (1995). 
R. Berkovits and Y. Avishai, Phys. Rev. Lett. 76, 291 (1996). 
D. Yoshioka and H. Kato, Physica B 212, 251 (1995). 



G. Bouzerar and D. Poilblanc, Phys. Rev. B 52, 10772 (1995); preprint cond- 
mat/9512153 . 



H.F. Cheung, K. Riedel and Y. Gefen, Phys. Rev. Lett. 62, 587 (1989) 
B. L. Altshuler, in Nano structures and Mesoscopic Phenomena, eds. W. P. Kirk and M. 
A. Read (Academic, San Diego, 1992); see also the article by Y. Imry in the same book. 
G. Vignale, Phys. Rev. B 50, 7668 (1994). 

P. Nozieres, Theory of interacting Fermi systems, (Benjamin, New York, 1964), p. 286. 
See, for example, J. A. Cullum and R. A. Willoughby, Lanczos Algorithms for Large 
Symmetric Eigenvalue Computations, (Birkhauser, Boston, 1985). 
[30] Our simulations were performed with the help of a parallel code on a cluster of up 
to 30 work- stations. The code was implemented with the help of software written by 
J. A. Holm, Master-Slave Implementation fur Parallel Virtual Machines, (Universitat 
Gottingen report, 1995, unpublished). Typically we averaged over 10 4 realizations of 
the disorder. 



12 



[31] A. Volker and P. Kopietz, preprint |cond-mat/ 9608027 . 

[32] D. Braun and G. Montambaux, Phys. Rev. B 52, 13906 (1995). 

[33] B. L. Altshuler and B. I. Shklovskii, Zh. Eksp. Teor. Fiz. 91, 220 (1986) [Sov. Phys. 

JETP 64, 127 (1986)]. 
[34] In Ref. 



11 



Altshuler, Gefen and Imry claim (see the text before their Eq.9) that for 
the calculation of £ 2 (/z, <£>) — S 2 (/x, 0) it is sufficient to consider the correlator K%(e, e', if) 
only for energies in the vicinity of the Fermi energy. For non-interacting electrons this 
is certainly not the case ||31||, so that one has to invoke inelastic processes in order to 



justify this assumption. 
[35] Ya. M. Blanter, Phys. Rev. B 54, 12807 (1996). 



13 



FIGURES 
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FIG. 1. Non-interacting average density of states of a 20 x 5 x 5 system with w = 4. For our 
numerical calculations we discarded all energies outside the interval [fi — Eq,/j] with /x = 0, Eq = 2 
(see the dashed lines) in order to obtain a model with approximately constant density of states. 
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FIG. 2. Probability distribution P(N) of the particle number in the simplified capacitance 
model on a 20 x 5 x 5 lattice for different values of e 2 /(CoA), with A -1 = po(fJ-) = 67. (a) 
e 2 /(C A) = 0; (b) e 2 /(C A) = A^/lOO = 0.67; (c) e 2 /(C A) = 3A-7IOO = 2.01; (d) 
e 2 /(CoA) = A _1 /20 = 3.35. Here and in all subsequent figures we have chosen w = 4, fi = 
and N + = 130.2. In this case N+ agrees with the average particle number, so that on average the 
ensemble of rings is not charged. 
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FIG. 3. Ratio of the even Fourier components p( 2m ) /P^ 2m ' (see Eq.^l]) of the particle number 
variance for different values of m as function of e 2 /(CoA). The solid line is the RPA prediction 
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FIG. 4. Zeroth Fourier component (E) of the variance a 2 (E, [/,, ip) as function of E/E c with 
E c = 1/180 and A" 1 = 67 for different values of e 2 /{C A). (a) e 2 /(C A) = 0; (b) e 2 /(C A) = 0.67; 
(c) e 2 /(C A) = 6.7. The band edge is at E /E c = 360. 
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FIG. 5. Second Fourier component p^ 2 \E) of the variance a 2 (E, fi,<f) as function of E/E c 
for different values of e 2 /(C A) with A" 1 = 67. (a) e 2 /{C A) = 0; (b) e 2 /{C A) = 0.67; (c) 
e 2 /(CoA) = 6.7. The solid line is the Altshuler-Shklovskii prediction ^5] with E c = 1/180 and 
f = 3.85. 
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